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ABSTRACT 

In this letter we examine the role of the first supernovae in proto-galaxies, their role in feedback and the 
consequences for disk assembly. Extending the picture proposed by Dekel & Silk (1986), we argue that 
energetic supernovae winds can expel baryons from all proto-galaxies with varying degrees of efficiency. 
The fraction of baryons retained and hence available to assemble into the baryonic disk is therefore, a 
function of the central velocity dispersion of the halo. Such a coupling of the baryonic component to the 
dark halo leads to the following interesting consequence, a prediction for a weak scaling of the zero-point 
of the Tully-Fisher relation or alternatively, the mass-to-light ratio with the central velocity dispersion of 
the halo. On applying to the case of the Milky Way halo, this feedback mechanism implies: (i) that the 
Milky Way halo lost ~ 10% of its original gas content, and the gas mass lost is roughly what is estimated 
for the mass in our X-ray halo consistent with the X-ray background in the soft band; (ii) a range in the 
inferred redshift of formation Zf, and the local baryon fraction fb for the Milky Way that depends on the 
initial spin parameter Xh of the halo. In a range of viable cold-dark matter cosmological models, we find 
that for a low spin halo (A/j ~ 0.02) - zt < 1, fb ~ 2%; for a median spin halo (A/i ~ 0.05) - Zf ~ 1 — 2.5, 
fb ~ 5%; and for a high spin halo (Ah, ~ 0.2) - zj ~ 4 — 8, fb ~ 20%. The observationally determined 
ages for the oldest disk stars in the Milky Way seem to rule out a low value for the spin parameter. 
Given the shape of the spin distribution of halos obtained in N-body simulations, while a high value of 
the spin parameter is not very probable, it is interesting to note that if this is indeed the case for the 
Milky Way halo, then feedback processes can cause the local baryon fraction to differ significantly from 
the universal value. 

Subject headings: galaxies — formation, fundamental parameters — galaxies: ISM — galaxies: 
kinematics and dynamics — galaxies: stellar content 



1. INTRODUCTION 

Despite tremendous recent progress both in the obser- 
vations of galaxies at high and low redshift and in the the- 
oretical understanding of their formation and evolution - 
there arc still gaps in our understanding of several of the 
physical processes that operate, most notorious of them 
being 'feedback'. The physics that couples baryons to the 
dominant dark matter component is comparatively poorly 
understood. There is compelling evidence from the many 
observed correlations in galaxy properties for the existence 
of feedback processes. In this work, we examine the feed- 
back provided by the first supernovae in proto-galaxies 
and the consequences for the subsequent disk assembly in 
the halo. In § 2 the feedback process is modelled, it is 
incorporated into a simple disk formation picture in §3, 
the consequences of the feedback are explored in §4 with 
emphasis on implications for the Milky Way disk and the 
conclusions are presented in § 5. 

2. FEEDBACK INDUCED BY EARLY SUPERNOVAE 

The very first SN to explode in a proto-galaxy have im- 
portant effects both in terms of setting into place feedback 
processes for subsequent episodes of star-formation as well 
as causing the initial expulsion of gas from the potential 
well. In this calculation, it is assumed that the first episode 
of star formation in a proto-galaxy occurs well prior to the 
assembly of the disk. This episode is expected to occur in 
small, isolated and dense clumps in the cold phase of the 
gas. 



The role of SN feedback and its implications for the evo- 
lution of a galactic nucleus have been discussed by Dekel & 
Silk (1986) [DS86 hereafter], Haehnelt & Rees (1993) and 
Haehnelt (1994). Haehnelt & Rees (1993) estimate the 
relative importance of several issues in the energetics: the 
impact on the heating-cooling balance due to the thermal 
energy injected by the SN; the inhibition of the collapse 
process produced by the radiation pressure from the ion- 
izing radiation as well as the mechanical energy of the ex- 
panding remnants; the effect on further star formation in 
the nucleus due to the consequent complex thermal struc- 
ture of the gas, and the production and mixing of metals. 
Semi-analytic models of galaxy formation (White & Frenk 
1991; Cole et al. 1994; Kauffmann, White & Guiderdoni 
1993; Kauffmann 1996; Baugh et al. 1997; Somerville & 
Primack 1998) have also successfully incorporated both 
the energy injection and mass loss due to supernovae. 

We focus here primarily on the dynamical effects of the 
first few SN. While the radiation induced ionization ef- 
fects act very locally, and on small scales ~ 100 pc, the 
energetic SN driven winds act on much larger scales since 
typical velocities range from ~ 5000 - 10,000 km s _1 . The 
mechanical effects of these winds are invoked as a means 
to expel baryons from the potential well prior to disk as- 
sembly. We show that this mechanism can modulate the 
fraction of baryons that are retained in a proto-galactic 
potential well and thereby naturally provides a variation 
in the fraction of available baryons to form a disk that is 
a function of the velocity dispersion of the halo. Addi- 
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tionally, SN explosions during these early epochs may ex- 
pel preferentially high angular momentum material and, 
therefore, the ejection could also aid in the subsequent ef- 
ficient disk assembly and the formation of central objects. 

First, we briefly outline the calculation in DS86 on the 
expulsion of gas from dwarf galaxies by supcrnovae driven 
winds. They computed the full energetics and the evolu- 
tion of the remnants through the two main phases, the adi- 
abatic expansion phase wherein the radiative cooling losses 
are estimated assuming radiative He + , Lyman-a and O2 
line-cooling, with a cooling function A(T) oc T^ 1 p 2 , for a 
metal-poor gas [Z <~ 1CP 2 Zm q ]', and the expansion into 
the momentum conserving snow-plow phase. Balancing 
the energy budget, they estimate a critical velocity disper- 
sion below which very little gas is retained in the poten- 
tial wells of dwarf-scale galaxy halos. This critical velocity 
was found to be v cr it < 100 km s -1 . Here, we take a more 
phenomenological approach to estimate the fraction of gas 
that can be expelled from halos with velocity dispersions 
exceeding this critical value. We assume that SN winds ac- 
celerate a variable fraction of the gas to the halo's escape 
velocity. Unless the initial mass function for this episode 
of star formation is considerably different, it is expected 
that SN would produce a patchy ISM leading to a complex 
relation between the fraction of baryons retained and the 
depth of the potential well. The crucial point is that even 
for massive halos a small fraction of the gas mass can be 
lost/expelled. We assume in accordance with DS86 that 
<~ 2% of the total SN energy output is available to drive 
the wind. The simple energetics argument used here is as 
follows: an upper bound on the fraction of gas mass lost 
can be computed by equating the escape velocity of the 
halo with the bulk velocity of the expelled gas. The total 
energy available per SN to drive the wind is given by: 

£ avail = 2 x 10 49 (_^_) erg. (1) 

V 0.02 M 10 51 erg 7 B V ' 

The total number of SN is A^n, and energy balance implies 
that, 
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where M; is the gas mass that is lost and w w i n d the wind 
velocity The escape velocity of the halo v csc — 2a; there- 
fore, for expulsion of the gas, we require that v^ ind > A a 2 . 
The total number of SN is proportional to all the mass that 
can potentially be turned into stars, i.e. the baryonic mass 
of the halo M^ ~ M^ , where is constant and is of 
the order of 5% in consonance with values determined from 
the standard picture of light element production from Big 
Bang Nucleosynthesis. Therefore, the total baryonic mass 
lost can be rewritten as, 

Mi = - oc - (3) 
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and on substituting the condition for gas expulsion, the 
fraction of baryonic mass lost /; is given by, 
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which can be translated into scaled units using the cali- 
bration obtained by DS86, 



/, = 5.0 x 10-' (-^-) ( 
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0.02 ' v 10 51 erg ' v 200 kms 
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By definition, it follows that the fraction of baryons re- 
tained f r = (1 — /(). Note that here we have assumed 
that the time taken for the first stars to go SN, expand 
and expel the gas via their winds is short compared to 
the time taken to assemble the disk. An attempt wherein 
feedback has been numerically implemented in a tree-SPH 
simulation of disk galaxy formation is reported in the re- 
cent work by Sommer-Larsen & Vedel (1998). Nulsen & 
Fabian (1997) have also examined the role of the SN feed- 
back in the context of the variation of the Faber- Jackson 
relation for low-mass halos. 

3. INCORPORATING FEEDBACK INTO SIMPLE DISK 
FORMATION MODELS 

Mo, Mao & White (1997) [MMW97] examined the for- 
mation of galactic disks in the context of cosmological 
models with a view to explain and understand the possible 
physical origin of their observed properties. Here we fol- 
low their notation and incorporate the above SN feedback 
into disk formation for 3 viable variants of the cold-dark 
matter (CDM) family of cosmological models - standard 
CDM [SCDM]: Q = 1, A = 0; Open CDM [OCDM]: 
O = 0.3, A = and ACDM with fi = 0.3, A = 0.7. 
The galactic-scale collapsed dark halo is modelled as an 
isothermal sphere truncated at an outer radius r2oo where 
the density is roughly 200 times the critical closure density 
Pcrit- A halo of a given mass M^ is fully specified by its 
velocity dispersion a and its redshift of collapse. The mass 
of the halo is, 



M h 
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where H (z) is the Hubble constant at a given epoch z, 
defined as, 

H(z) = 100ft [fi A + (1 - fiA - flo)(l + zf 

+ fi (l + zf] 1 / 2 kms- 1 Mpc- 1 (7) 

using the same notation as MMW97. The total bary- 
onic mass available to form the disk, Md; s k oc f r fi, M na io, 
therefore from eqns. 6 and 7, 

M disk - 2.63 x 10 11 ft" 1 [1 - 5.0 x lO" 2 (— f )- 2 ] 

200 kms 1 
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Initially the baryons and the dark matter have comparable 
torques but dissipation causes differential settling of the 
baryons. We assume that the total angular momentum of 
the disk J d is proportional to the angular momentum of 
the halo Jh, Jd = fj Jh- The angular momentum of the 
disk is computed assuming a flat rotation curve neglect- 
ing the inner-most region where the disk is self-gravitating 
(models that take into account the back-reaction of the 
halo have been worked out in detail by Blumenthal et al. 
(1986), Dalcanton et al. (1997) and MMW97). The spin 
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parameter A/j of the halo is a measure of the acquired angu- 
lar momentum during the collapse process, and is defined 
as \ h = J G^ 1 M~i where J, E and M are the total 
angular momentum, energy and mass of the halo. Assum- 
ing that the disk settles down to an exponential radial 
surface density profile, the following expression can be ob- 
tained for the disk-scale length r d in terms of the relevant 
parameters, 



A 



fi 



r d = 8.17/ l - 1 (^)(^)[l-5.0xl0- 2 ( 

x(^)-kpc. ^ 
no 



200 kms 
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(9) 



The Tully-Fisher (T-F) relation can be written as, 
log [L disk ] = log [L ] + (3 log 



200 kms- 



(10) 



where (3 is the Tully-Fisher index and L the zero-point 
of the relation. A recent determination by Giovanclli 
et al. (1997) in the I-band, yields (3 — 3 and Lq = 
3.5 x 10 10 h~ 2 Lq. Although the value of the index (3 does 
depend on the wave-band {(3 = 2.5 in the B-band), the 
scatter is small (roughly 0.35 mags.) is magnitude depen- 
dent and increases steeply with decreasing velocity width. 
Assuming an efficiency ratio T for the conversion of gas 
in the disk into stars, the disk mass Md (eqn. 8) can be 
translated directly into a luminosity. Comparing this with 
the T-F relation, we require, 
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There are now 2 possible ways to interpret the above equa- 
tion, in terms of (i) yielding a zero-point Lq for the T-F 
relation that does depend on the depth of the potential or 
(ii) if we require Lq to be independent of a, then it follows 
that T is a weak function of the velocity dispersion of the 
halo. Even if the T-F zero-point depends on a, we find 
that it is independent of the spin of the halo Ah as first 
pointed out by White (1987). 

4. CONSEQUENCES OF THE FEEDBACK 

A stronger coupling between the velocity dispersion and 
disk mass arises as a result of the feedback process. We 
make no claims with regard to understanding why the T- 
F exponent is 3, it is assumed to be so motivated by the 
observations. Any successful theory of galaxy formation 
needs to account for both the slope of the relation as well 
as the low scatter. 

4.1. Application to the Milky Way disk 

In the context of this feedback process, we can now ex- 
amine what the observed parameters of the disk of our 
Galaxy imply for its formation history. The age of the 
disk stars in the Milky Way can be used as a diagnostic to 
pin down the redshift of assembly. Also, constraints can be 
obtained on the frequency of accretion events in our halo 



if the bulk of the stars in the Milky Way disk assembled 
at z > 1. The thick disk stars have a relaxed exponential 
profile and were hence probably in place quite early on. 
Moreover, since the relaxation time-scale for stellar sys- 
tems is long compared to their typical dynamical times, 
these stars most probably formed in situ. Chemical evo- 
lution models of our galaxy seem to indicate that a con- 
stant star-formation rate is not a bad description (Wyse 
& Gilmore (1995)), implying that half the stellar mass of 
our Galaxy probably assembled well before z = 1. 

An early generation of stars (progenitors of these SN) 
also provides a natural solution to the so-called G-dwarf 
problem, i.e. the paucity of metal poor stars in our 
Galaxy that closed-box chemical evolution models predict 
in abundance (Pagel 1989). Observationally determined 
parameters for our Galaxy are: a disk mass AIj\/[w ^ 
6 x 1O 1O M , velocity dispersion a <~ 155 km s _1 and 
a disk-scale length <~ 3.5 kpc. Firstly, in accordance 
with this prescription, the predicted gas fraction lost by 
the Milky Way halo, fi ~ 10%. Substituting the appro- 
priate values in eqns. 8 and 9, assuming maximal efficiency 
for disk formation = fj and H = 50 km s _1 Mpc -1 , 
we obtain the following relations: 



/ x h , H{z) ! 
l 0.05 M H ' 
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Combining these two bounds yields, 



(_^L)(JL)-i = 1.054. 
V 0.05 M 0.05 ; 



0.26(712) 



(13) 



The constraint above can be interpreted as providing 
bounds on the value of the local baryon fraction depend- 
ing on the value of the initial spin parameter of the Milky 
Way halo. Numerical simulations and analytic estimates 
seem to indicate that the spin distribution of halos P(X) 
can be well-fit by a lognormal with a mean value of 0.05 
(Ryden 1988; Warren et al. 1992). If the spin of the Milky 
Way halo was low A^ <~ 0.02 then eqn. (13) implies that 
the local baryon fraction (/& ~ 2%) is lower by a factor 
of two compared to the universal value; a mean value for 
the spin parameter yields fb ~ 5% in agreement with the 
universal value whereas a high value of the spin parameter 
(A/j ~ 0.2) implies /& <~ 20% which is much higher than 
the universal value but is in good agreement with the value 
inferred from X-ray studies of the Intra-cluster medium in 
clusters of galaxies. 

A feedback scenario such as the one suggested in this 
work - in which only 10% of the gas is ejected has interest- 
ing implications for the origin of the extended X-ray halo 
around our Galaxy. The limit on the gas mass expelled 
from the halo of our Galaxy predicted by this model is in 
good agreement with the mass of the hot corona of our 
Galaxy (Bregman 1980; Kerp et al. 1998) consistent with 
the soft X-ray Background. Given values of fb and Xh eqn. 
(12) also then provides naturally an estimate for the for- 
mation redshift for our Galaxy - (i) for a low spin, /(, ~ 2% 
Zf < 1, (ii) a mean spin, fb ~ 2% Zf ~ [1.0 — 2.5], and 
(iii) for a large spin, f b ~ 20% z f ~ [4.0 - 8.0], for the 
3 cosmological models considered here [as shown in Fig. 
1] . The ages of the oldest disk stars in our Galaxy (Wyse 
1997), seem to be consistent with a high- redshift of for- 
mation, thereby ruling out a low spin parameter for the 
Milky Way halo. 
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5. DISCUSSION AND CONCLUSIONS 

In this letter we have examined the consequences of 
feedback from an early episode of SN explosions in proto- 
galaxies. An extension of the mechanism proposed by 
Dekel & Silk (1986), to higher velocity dispersion halos 
yields a modulation of the baryon fraction available for 
disk formation that depends on the depth of the poten- 
tial well. This feedback applied to simple prescriptions 
for disk formation implies that the either the zero-point 
of the Tully-Fisher or the mass-to- light ratio is weakly 
dependent on the velocity dispersion of the halo. The im- 
plications of this scenario for the Milky Way disk are that 



approximately 10% of the original gas that fell in was ex- 
pelled prior to disk assembly, a high-redshift of formation 
which is in good agreement with the ages of the oldest 
thick disk stars is inferred for average to high values of the 
spin parameter of the Milky Way halo with a local baryon 
fraction that ranges from the universal value of 5% to a 
value of 20% which is what is inferred from X-ray studies 
of the Intra-cluster medium in clusters of galaxies. 

Useful conversations with Martin Haehnelt, Martin 
Rees, Steinn Sigurdsson, Joe Silk, Rachel Somerville and 
Rosie Wyse are gratefully acknowledged. PN thanks the 
anonymous referee for very constructive comments. 
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Fig. 1. — The constraint on the rcdshift of assembly for the Milky Way disk obtained from early SN feedback is marked by the hatched 
region. The solid curve is the SCDM model, the dashed curve OCDM and the dotted curve is for the ACDM model. 



